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Abstract. 

The random vector potential model describes massless fermions coupled to a quenched ran- 
dom gauge field. We study its abelian and non-abelian versions. The abelian version can 
be completely solved using bosonization. We analyse the non-abelian model using its su- 
persymmetric formulation and show, by a perturbative renormalisation group computation, 
that it is asymptotically free at large distances. We also show that all the quenched chiral 
current correlation functions can be computed exactly, without using the replica trick or the 
super symmetric formulation, but using an exact expression for the effective action for any 
sample of the random gauge field. These chiral correlation functions are purely algebraic. 
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Two dimensional random systems are much less understood than their pure companions. 
It is not clearly understood how (if possible) conformal field theory techniques can be applied 
to random systems. For example, what kind of (probably non-unitary) conformal field 
theories could describe the infrared fixed points of two dimensional random systems is still 
an open question. 

Besides the familiar study of the effects of disorder on phase transitions [l], |], |[ |J , there 
has been a renewed interest in two dimensional random systems in connection with the 
quantum Hall phase transition ||, |6|, |7|, or in connection with models of randomly pinned 
flux lines in superconductors f§. 

Obviously, we will not give a complete answer to the question formulated above. However, 
the models we will study are simple enough, and possess enough symmetries, to allow for 
exact computations which do not rely on the replica trick or on the supersymmetric formu- 
lation. Using this direct computation, we derive some informations on the operator product 
algebra governing the infrared behavior of these random models. Although our result are 
only partial results, we hope they provide an information useful enough to start answering 
the question. 

The abelian version of the models was recently introduced for analysing the quantum 
Hall transition in ref. [01. The non-abelian version is similar to the model introduced in ITU 



for describing the effects of impurities interacting with fermions close to the Fermi surface 
in 2 + 1 dimensional systems. However, as explained in the Appendix, it seems that our 
conclusions are in disagreement with some of those of ref. 



1 The model. 



i) Definition of the model. 

The model describes N massless Dirac fermions minimally coupled to a random non- 
abelian gauge field with euclidean action : S = - J d 2 x^(i0 +A^)ty. It is useful to introduce 

the complex coordinates z = x + iy,~z = x — iy and the components of the fermions (?/>+, ip 3 + ) 
and J = A^. The action then becomes : 



S[A] = J ^ 0- ; ; (dzSl + AQ ^ + (dJt + A{ 



—k 



where A\ k = iJ2 a A°;(t a y k , with (A a z )* = A|, is the gauge field. Here the hermitian matrices 
t a form the iV-dimensional representation of U(N). We denote by f abc the U(N) structure 
constants : [t a ,t b ] = if abc t c . The Dirac fermions take values in this iV-dimensional represen- 
tation. 

The gauge field is assumed to be a quenched variable with a Gaussian measure : 



P[A] = exp 



a J n 



(2) 



We will be interested in computing the quenched average of the correlation functions of 
the currents J® = (\E'7 At t a \l/). Explicitely, their components are : 



ja 



(3) 
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A convenient way of doing it consists in introducing sources for the currents. This amounts 
to shift the random field A by an external sources a source : A — > A = A + H SO urce- We are 
then interested in computing the average of the logarithm of the partition functions with 
sources : 



T[a SOUIce } = w[A] = h g z[A] (4) 

with A = A + Hsource and 

Z[A] = e w[A] = J exp (S[A}) = Det[0 +4]. (5) 
W[A] is the generating function for the connected Green function : 



!«(y)-)J 



5W[A) 



6A b u {y)SA a Jx) 



ii) The abelian case. 

Before plunging into the non-abelian case, let us first present a simple way of solving 
the abelian case which corresponds to AT = 1. This model was introduced in ref.0 in 
connection with the quantum Hall transition, and studied there using the replica trick. For 
N — 1, the action (Jl]) can be bosonized. Using the standard rules, = e^d^, and 

= |(<9 M $) 2 , one finds the bosonic form of the action ([!]) : 

d 2 x (\ 



S N=1 [^,A} =J^ (l^f + iA^d^ . 



Since it is a gaussian model, it can be easily solved without using the replica trick. Let us 
introduce the Hodge decomposition of A^ : A^ = <9 M £ + e^dyf]. The fields £ and 77 decouple 
in the measure (Q) : 



P[£;rj\ = exp 



d 2 x 



71 



(6) 



The action 5' iV=1 [$,y4] is independent of £ and therefore the field £ is irrelevant. This 
fact was expected since the field £ represents a pure gauge whereas the physically relevent 
quantity is the field strengh F = e^ v d^A u = (d^d^rj. Moreover, the field r\ can be absorbed 
into a translation of $ : 

S N=1 [$, A, = e, u d vV ] = S N=1 [<S> + t7] ,A = 0] + y — (d, V ) 2 . (7) 



7T 



This equation encodes the anomalous transformation law of the determinant Det[i0 
under a chiral gauge transformation of the abelian gauge field A : 



Detlih 



Det[0\ 



exp 



The fact that the field rj can be absorbed into a shift of $ does not mean that the 
quenched correlation functions are identical to those in the pure system. Using eq.(^), we 
have 



(II e^^U^d* = aM (II e-™*^)) A=0 . 
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It factorizes into the product of the correlation functions in the pure system times a simple 
function of the impurities. However, the average of this function is not irrelevant since the 
variables rj have long-range correlations : r)(x)r](y) = — ncr^d^d^)' 1 ^, y) = — cr log \x — y\. In 
particular, it changes the values of the critical exponents. The dimensions A of the vertex 
operators exp(ia$(x)) in the quenched theory are : 

a ct _ a a _ ^-^2 

quenched pure ^ " 

Averages of product of correlation functions can be computed similarly. One has : 



n m 

x (Y[e ianHxn) )A=o ■■■ (n e ^* (ym) >A=o 

n m 

where ((•■■)) refers to the //-correlation functions with the free field measure @. Clearly, 
conformal invariance is unbroken in the random abelian case. More details will be described 



elsewhere [12 



This way of solving this very simple model is closely related to the change of variables used 



in ref.|TT|. What we are going to present is the non-abelian analogue of this construction. 



In the following, we will restrict ourself to the SU(N) sector of U(N). 

Hi) Symmetries and Ward identities. 

The currents J® = (\E'7 /J t a \l/) are conserved in the pure system, i.e {(d^J^) 0)a=o = 0. 
It is well known in conformal field theory that in the pure system these currents represent 
a Kac- Moody algebra, alias a current algebra, of level one This is non longer true in 
the presence of impurities. However, the measure -P[A], eq.(|2|), is invariant under a global 
SU(N) rotation : A a — > UA^U^ 1 with U G SU(N). The Noether current associated to this 
symmetry is A u , which is therefore conserved inside any quenched correlation functions. I.e. 
{d^Afj) (Oi)a- ■ ■ (Om)a — 0. By integration by part, this is equivalent to the conservation 
law for the currents J£ defined by : 
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V = • (8) 



8A* 



Explicitely, the Ward identities are : 



d, {J?[(O 1 ) a ---{OmM)=0- (9) 
Equivalently, by computing the action of J* on the expectation values, we get : 



E (°i)a ■ ■ ■ [(fl^i - (d^) A (o 3 ) A \ ■ ■ ■ (o M ) A = o. 

j 

Therefore, in the quenched theory the conserved currents are not but their insertions in the 
connected Green functions. This point can also be recovered using the replica formalism. 
(There also exist Ward identites for chiral gauge transformations : A z — > UA z U~ l and 
A Y — > VA-V^ 1 ). As usual, the Ward identites (|9[) are up to contact terms which encode the 
transformation law of the fields. 

Usually, conformal field theories with a global SU(N) symmetry group possess two chiral 
sets of SU (N) conserved currents satisfying a Kac- Moody algebra. So, it seems reasonable 
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to expect that the field theory describing the large distance behavior of this quenched model 
should possess two sets of currents. We will verify this point in the following section while 
computing the chiral quenched correlation function. We will actually obtain a little more 
since we will find that the chiral correlation functions are purely algebraic and do not present 
any crossover between their ultraviolet and infrared behaviors. 

2 The supersymmetric approach. 

A standard approach for studying this random model would be to use the replica trick or 
the supersymmetric formalism. In this section, we use a perturbative renormalization group 
computation to derive the infrared behavior of a few averaged correlation functions. 

i) The supersymmetric action. 

Since the partition function is a determinant, we choose to use the supersymmetric formal- 
ism. Therefore, we introduce auxiliary bosonic fields, denoted (6j,c J ), in order to represent 
the inverse of the partition function ([5]) : 




DbDc exp(-S aux [A}) 



with 

S aux [A] =j*± (b, (ckSi + AQ c k + b 3 (d z 5{ + A{ k ) c k ) (10) 

In order to calculate averages of products of correlation functions, we need to introduce as 
many copies of the fermions and the auxiliary bosons as necessary. We denote them by 
(^ a , and c a ). Their dynamic is governed by the action St^LA] with : 

S tot [A] = J2(S a [A] + S: ux [A}). (11) 

a 

The fields are only coupled through the impurities. After integration over A with the measure 
(0) we obtain the effective action : 

Seff = E (S a [A = 0] + S: UX [A = 0])-a J £ f (£#T)- (12) 

a a a a 

where H% and H§ denote the total SU (N) currents : 

K = ^-Ani^ + b^ty^ 

The action S e ff is explicitely supersymmetric. More precisely, the action St f[A] is su- 
persymmetric before any integration over A. There are two supercharges, Q and Q, acting 
separately on the left and right sectors. The charges Q, acting only on the left movers, is 
defined by : 

Q(V-) = ; Q(b) = 
Q(V+) = c ; Q(c) = 0. 
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Q acts trivially on A and on the right sector. Clearly, it satisfies Q 2 = 0. One verifies 
that S'totfA] is susy invariant : Q(S to t[A]) = Q(S tot [A]) = 0. One also verifies that the total 
current H z and H§ are respectively Q and Q closed. This means that there exist two local 
fields and K% such that : 

H a z =Q(K a z ) and H§ = Q(K§). 

This property ensures that the total partition function, obtained by integrating over all the 
fields (the fermions plus the bosons), is independent of A as it should be. 

There exists an intriguing similarity between the above construction and conformal topo- 
logical field theory in two dimensions. In the latter, one uses a twisted version of a N = 2 
supersymmetry to construct the theory such that the stress tensor is susy closed. The closure 
of the stress tensor ensures that the partition function of the topological theory is indepen- 
dant of the background metric. In the random model, we can extend the supersymmetric 
algebra we just described such that the new algebra possesses two supersymmetric charges. 
The fact that the current is susy closed is the analogue of the fact that the stress tensor is 
susy closed. It seems natural to wonder if techniques of conformal topological field theories 
cannot be translated into new methods of analysis of two dimensional random systems. 

ii) A renormalisation group computation. 

The only coupling constant which can be renormalized is a. As usual, the one-loop j3 and 
7-functions are encoded in the operator product expansions of the fields in the unperturbed 
theory. These can be read off from the action S to t[A = 0]. We have : 

■0+(z) ip--j(w) ~ - — - — r + regular, 



[z — w) 
5' 



c (z) bJw) ~ - — - — r + reqular. 
v ' n ' (z-w) 

The total currents H z satisfy a Kac-Moody algebra with zero central charge : 

• .cabc 

H a z {z) H b z (w) ~ -i— H c z {w) +reg. 

There is an exact compensation between the central charges associated to the fermions and 
to the auxiliary bosons. 

The currents H z have dimension one; therefore, the f3 function vanishes at tree level. At 
one loop, one finds : 

d = /3(a) = -C G a 2 + CV), (13) 

where Cq > is the Casimir in the adjoint representation, f abc f dbc = CG$ ad . The sign in 
eq. ([l3|) is important; it tells us that the model is asymptotically free in the infrared regime. 
We have : 

/TV) = 

a[ ] 1 + a C G ]og(R/a) ' 

with cr > the value of the coupling constant at the lattice cut-off a. Since the model is 
asymptotically free in the infrared, we can use the renormalization group to evaluate the 
large distance behavior of the two-point correlation functions using the usual formula : 

(HRMO)U = (Ha)*(o)Um ™p (-2 dw jj=y) > ( 14 ) 
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where 7$ is the 7-function for $. 

Let us first compute the asymptotic of the chiral two-point function (J%(z) J^(0))a- To 
evaluate it we only need one copy of the fermion and their supersymmetric partner. We find 
that the 7-function of J" is 7j(cx) = 1 + 0(a 2 ); i.e. there is no one-loop correction even 
though are not conserved. Therefore, at large distances we have : 

(15) 

In particular, there are no logarithmic correction. 

Consider now the average of the product of two such two-point functions but of oppo- 
site chiralities, i.e. {J^(z)J^(0))a{J§{z)J^(0))a- To compute it we need two copies of the 
fermions. These quenched averages are then represented by the two-point functions of the 
operators = J^ l \z)J^ {z), products of currents in the different copies. For these op- 
erators the 7-functions have one-loop corrections. There is a mixing between the operators 
with different indices : 

= 26 a b 6l + 2aftff + 0(a 2 ). 

In particular the SU(N) scalar operator O = J2 a J^'^Jj does not mix with the others, 
and 7 = 2 + 2C G + O (a 2 ). Therefore, 



E (j«(z)jm)A(j§(z)jm)A ~ ^ (log \zz\y m 

0,6 \ ZZ 7 

Notice the logarithmic correction which breaks the self-averaging property. It is remarkable 
that the logarithms appear in the averages of products of correlations of opposite chiralities, 



eq . ([lq) , but not in the averages of the chiral correlations, eq . (|15|) . 

Infrared asymptotics of averages of higher products of current two-point functions, or of 
fermion correlations, can be computed similarly. In particular, we find logarithmic correc- 
tions for the quenched two-point function (e(i?)e(0))A for the energy operator e = : 

,2 \ 



(e(z)e(0)) A ~ - (log \zz\) 



-±c v /c G 



zz 



where Cy is the SU(N) Casimir in its iV-dimensional representation. 



3 The quenched chiral correlation functions. 

In this section, we describe how some of the quenched correlation functions can be directly 
computed without using the replica trick or the supersymmetric formulation. This relies on 



the Polyakov-Wiegman (PW) formula |14| for the effective action W[A], eq.flj). Once this 



formula has been recalled, the computations are rather simple. 

i) The Polyakov-Wiegman formula. 

The PW formula relies on the fact that in two dimensions the effective action can be 
exactly computed by integrating its anomalous transformation under a chiral gauge trans- 
formation. This is the non-abelian analogue of eq.(|7]). Indeed, let G M = Yjat a G a AA\ be the 
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generating functions of the connected current Green function in the pure system : 

G a M\ es {j;) A = vr^l. (17) 



They satisfy the anomalous Ward identities HT4f , 

d z G-+chG z + [A z ,G-} + [A-,G z } = 0, 

d z G?-dzG z + [A z ,G?]-[A- z ,G z ] = 2F ZZ [A], (18) 
with F ZZ [A] = d z A- - c\A z + [A g , A-}. Eqs.([T|) completely specify G^A] : 

G ' [A] = A >-lk^A^ (19) 

G ^ = A ^8^ZT^ A ' (20) 

Here ad.A z denoted the adjoint action of A z . Notice that G 2 [v4] is local in A z but non-local 
in Ag. 

ii) Explicit formula for the chiral correlation functions. 

Knowing explicitely the correlation functions for any impurity sample, it is a priori pos- 
sible to take the quenched average. But let us first concentrate on the average of the chiral 
correlation functions involving only currents of the same chirality; e.g. involving only J". 
Consider first the average of products of one-point functions ([H])- Since the quenched av- 
erage is defined by A®(x)A^(y) = airS ab S^(x — y) and since (J z )a are linear in A z , these 
quenched correlations can be computed using the Wick theorem applied on A. For exemple, 
the two-point and three-point functions are : 



(J«{z,)) A {J b {z 2 )) A = 2a*5 ab {j_9 z ^ 

(j«( Zl )} A (J b z (z 2 )} A (J c M))A = i{™yr bc (^d z ) 



lob' 



nb 



3a 



+ (cyclic permutation) 

2 

(z x - z 2 )(z 1 - z 3 )(z 2 - z 3 ) 



a l \ a l 

L/ z / (zz,z\) VL^z/ («3,«a). 



More generally, the average of products of one-point functions is the sum of connected 
correlations which can be expressed in terms of the correlation functions of the pure system : 



(JHzi))a ■ ■ ■ (J^(z M )) A \ C ° nneCted = Mc M - l (J?{ Zl ) ■ ■ ■ J a z M (z M )) (21) 

Here, (■ • -)o denote the pure correlation functions. They are known exactly ]TJ|. The relation 
between the quenched correlations and their connected parts is the usual one. 

More interesting is the average of products of correlations with insertion of the conserved 
currents (§) since they encode the underlying symmetry algebra. We find : 

J^( Zl ) ■ ■ ■ J?"{z M )(J?{wi)J${£l) ■ ■ -ft ■ ■ ■ (j;'(w P )J b s(tp) ■ ■ -) A 

= (j?( Zl ) ■ ■ ■ J n z M {z M )J a Awi)J h At,i) ■■■)*■■■ {JT(w P )j b z -a P ) ■ ■ -) 
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M 

+ E<^i) • • • ^ J (%) • • • Jz M M J? M • • -)o x (22) 

x (J?( Zj ) J?{w 2 )jf{i 2 ) • • -) • • • (J"' " " ')o 

+ ■■• 

+ (j? k) j, 6i (6) • • -> • • • (^(^i) • • • r(zM)jrK)^^) • • % 

The hatted fields have to be omitted. Here, we assumed that there is no insertion of one-point 
functions. The formula ( p2[) is actually simpler in words : it is obtained by distributing the 
currents J™ 1 ^), • • • , J™ M (z M ) among the pure correlators in all possible way, each counted 
only once. 

Notice that all the chiral quenched correlation functions are purely algebraic, without 
any logarithmic correction, in agreement with the renormalization group computation. So, 
conformal invariance does not seem to be broken in the random model. 

From eq.(P2"D we read the operator product expansion of the fields. The currents satisfy : 

J?(zi)J?{z*) = V2 + 1 J ns (z 2 ) + reg. (23) 

[zi - z 2 y zi - z 2 

Therefore, the quenched conserved currents satisfy the commutation relations of a Kac- 
Moody algebra, exactly as the currents in the pure system do. 

The operator product expansion between the conserved currents and the correlators 
{J^{w 1 )J^{w 2 )---) A are : 

+ 

These operator product expansions are unusual in conformal field theory with Kac-Moody 
symmetry. In particular, they imply that the fields {J^{w\) J^ 2 {w 2 ) ■ ■ -) C A are not associated 
to highest weight vector representations. We don't know to which category of representations 
they correspond to. 

Contrary to the chiral quenched correlation functions which are easy to compute, the av- 
erages of correlation functions involving fields of opposite chiralities are difficult to evaluate. 
This is due to the fact that the generating functions G Z [A] and are non-local in A- 

and A z respectively. A naive perturbative expansion is spoilled by untractable divergences. 
This could have been anticipated in view of the logarithmic corrections present in eq . (|16|) . 
The occurence of these logarithms probably means that the gluing of the left and right sec- 
tors is more subtle than usual. In particular, one has to learn how to deal with an infinite 
number of primary fields. We don't know any good algebraic way of doing it correctly in the 
non-abelian case, but we hope that the remarks presented in this note will be useful to find 
the answer. 



jf na i a ' 3 

E V? W • • • ^ K ) -~)a + ^9. (24) 

y z-w 3 
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4 Appendix. 



In this appendix, we explain why our results seem to be in contradiction with some of those 
presented in ref.jnj. The approach used in this paper is based on the replica trick. Thus, 
the authors introduce n copies of the N fermions, forming alltogether Nn Dirac fermions. 
These fermions form a representation of level one of the SU(Nn) Kac-Moody algebra. After 
averaging the random potential, these Nn free massless Dirac fermions are coupled by a 
current- current interaction : 

Sm = v [^(EJF)(£J?) (25) 

J 71 p p 

where the currents J^ ,p are the SU(N) currents in the p th replica. Not all the degrees of 
freedom are coupled by this interaction since only the diagonal sum of the SU(N) currents 
interacts. The authors of ref.||10|| argue that the interaction opens a gap and that some of 
the modes become massive. More precisely, let us decompose the SU(Nn)k=i representation 
as SU(Nn)k=i = SU(N) k=n x SU(n) k=N . They argue that the modes corresponding to 
SU(N)k =n become massive and therefore that the infrared behavior is described by the 
SU(n)k = N conformal field theory. In particular, this would implies that the correlation 
functions of the diagonal SU(N) currents (J2 P J^' p ) decrease exponentially at infinity. But, 
since these correlation functions at n = give the quenched averages of the connected 
correlation functions of the currents, the latter statement is in contradiction with the result 
we obtained using the Polyakov-Wiegmann formula. 
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